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Ecologists  have  long  debated  whether  insect-herbivore  population  dynamics  are 
driven  by  “top-down”  forces1,  such  as  predators  and  pathogens,  or  “bottom-up” 
forces2,  such  as  induced  plant  defenses.  In  the  case  of  outbreaking  forest  insects, 
theoretical  ecologists  have  argued  that  top-down  forces  are  most  important3,  be¬ 
cause  mathematical  models  that  incorporate  only  predators  and  pathogens  ex¬ 
plain  time  series  of  outbreaks  from  the  field4.  Proponents  of  bottom-up  forces 
concede  that  induced  plant  defenses  often  have  only  weak  direct  effects,  but  nev¬ 
ertheless  argue  that  the  models  that  do  not  consider  plant  quality  are  inaccurate 
because  they  do  not  include  interactions  between  induced  defenses  and  insect 
pathogens5,6.  Such  interactions  reduce  average  infection  risk  in  the  laboratory7, 
and  should  therefore  reduce  infection  rates  during  outbreaks  in  nature,  when  in¬ 
duced  defenses  are  highest8,  but  data  from  the  field  have  shown  that  infection 
rates  instead  peak  during  outbreaks9  11 .  Here  we  reconcile  models,  experiments, 
and  field  data,  first  by  using  a  field  experiment  to  show  that,  in  the  outbreaking 
North  American  gypsy  moth,  induced  plant  defenses  affect  pathogen  transmis¬ 
sion  mainly  by  reducing  variability  in  infection  risk  rather  than  average  infection 
risk,  which  has  the  effect  of  increasing  infection  rates  during  outbreaks,  not  re¬ 
ducing  them12.  We  then  show  that  this  reduction  in  variability  has  important 
implications  for  insect  population  dynamics,  because  models  in  which  variability 
in  pathogen  infection  risk  is  determined  by  induced  plant  defenses  produce  realis¬ 
tic  outbreaks,  whereas  standard  models,  when  accurately  parameterized,  instead 
predict  unrealistic  stable  population  dynamics.  Our  work  makes  clear  that  top- 
down  and  bottom-up  forces  may  interact  in  complex  ways,  contrary  to  traditional 
either-or  viewpoints,  and  suggests  that  induction  of  plant  defenses  may  augment 
the  effectiveness  of  baculoviruses  in  reducing  the  damage  caused  by  forest  pests13. 

Pathogens  of  outbreaking  insects  are  often  baculoviruses,  fatal  diseases  transmitted  when 
host  larvae  accidentally  consume  foliage  contaminated  by  virus  from  infectious  cadavers  of 
other  larvae14.  Infection  occurs  as  a  result  of  virus  consumption,  and  so  plant  defensive 
compounds  consumed  together  with  infectious  viral  particles  can  alter  the  infection  process 
within  the  insect6.  Consumption  of  plant  defensive  compounds  has  therefore  been  shown 
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to  reduce  average  infection  risk  in  the  laboratory6.  Because  defensive  compounds  are  often 
induced  by  defoliation15,  and  because  defoliation  and  thus  induced-defense  concentrations 
increase  with  insect  densities8,  the  laboratory  data  imply  that  bacnlovirns  infection  rates 
in  the  held  should  decline  with  increasing  insect  density7.  In  nature,  however,  infection 
rates  instead  increase  with  density,  causing  population  collapses  at  outbreak  peaks9-11.  One 
explanation  for  this  contradiction  is  that  conditions  in  laboratory  experiments  differ  greatly 
from  those  in  nature6,  and  laboratory  experiments  therefore  may  be  irrelevant  to  processes 
in  nature16.  To  date,  however,  there  have  been  no  convincing  experimental  tests  of  whether 
induced  defenses  affect  baculovirus  transmission  in  the  held. 

To  carry  out  such  a  test,  we  used  the  gypsy  moth  ( Lymantria  dispar ),  its  baculovirus,  and 
one  of  the  gypsy  moth’s  main  host  trees  in  North  America,  the  red  oak,  Quercus  rubra1' .  A 
previous  study  relied  on  experimental  defoliation,  without  successfully  causing  induction18. 
It  therefore  appears  that  defoliation  must  be  quite  severe  for  induction  to  occur,  yet  se¬ 
vere  defoliation  would  remove  so  much  leaf  material  that  it  would  be  impossible  to  measure 
virus  transmission  in  the  held.  Accordingly,  we  instead  induced  red  oak  defenses  by  spraying 
branches  with  jasmonic  acid  or  “JA”  19 .  JA  is  a  plant-signalling  compound  that  in  seeching 
red  oaks  increases  hydrolyzable  tannin  concentrations20,  an  induced  defense8  that  strongly 
affects  average  gypsy-moth  infection  risk  in  the  laboratory21.  Reassuringly,  our  JA  treatment 
induced  hydrolyzable  tannins  to  the  same  extent  as  defoliation  in  nature8  (Table  1,  note  that 
control  branches  were  sprayed  with  an  identical  solution  that  lacked  JA,  and  that  induction 
affected  only  sprayed  branches,  not  entire  trees,  see  Supplemental  Information).  We  then 
mimicked  natural  transmission  by  first  adding  virus-infected  larvae  to  the  branches,  allowing 
the  larvae  to  die,  and  then  allowing  uninfected  larvae  to  feed  on  the  branches  for  one  week, 
a  period  of  time  short  enough  to  ensure  only  one  round  of  transmission  (see  Supplemental 
Information  for  details). 

This  experiment  showed  that  the  main  effect  of  induction  is  that  it  reduces  variability  in 
infection  risk.  To  see  this,  we  used  the  mathematical  theory  of  epidemics12  to  quantify  the 
effects  of  induction  on  average  infection  risk  and  variability  in  infection  risk  in  onr  experiment. 
According  to  the  theory,  the  fraction  infected  i  in  onr  experiments  can  be  expressed  as, 

=  T  l°g  (l  +  C2vPaT)  . 
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log(l  -  i) 


(1) 


Here  T  is  the  length  of  the  experiment,  and  Po  is  the  virus  density,  which  is  constant  in  our 
experiments.  By  fitting  equation  (1)  to  our  data,  we  estimated  u,  the  average  infection  risk, 
meaning  the  average  instantaneous  transmission  rate,  and  (7,  the  variability  in  the  infection 
risk,  meaning  the  coefficient  of  variation  of  the  distribution  of  transmission  rates.  Variability 
in  infection  risk  is  thus  expressed  as  a  distribution  of  transmission  rates,  such  that  some  hosts 
have  high  risk  and  others  have  low  risk,  following  a  distribution  with  mean  v  and  coefficient 
of  variation  C.  To  see  the  effects  of  induction  on  variability  in  infection  risk,  note  that  when 
variability  C  >  0,  equation  (1)  predicts  that  —  log(l  —  i)  is  a  nonlinear  function  of  virus 
density,  whereas  when  C  — »  0  we  have  —  log(l  —  i)  —  vPqT,  so  that  —  log(l  —  i)  is  a  linear 
function  of  virus  density. 

The  reduction  in  variability  in  infection  risk  due  to  induction  is  then  visually  apparent  in 
our  data,  in  that  transmission  on  non-induced  branches  is  a  nonlinear  function  of  virus  density 
(fig.  1A),  while  transmission  on  induced  branches  is  linear  (fig.  IB).  Induction  thus  reduced 
variability  in  infection  risk  C  from  a  high  level  to  near  zero  (fig.  1C),  but  the  reduction  in 
average  infection  risk  was  relatively  small.  Average  infection  risk  and  variability  in  infection 
risk  were  similarly  lower  on  induced  foliage  in  the  laboratory  (Supplementary  Information), 
suggesting  that  the  effects  seen  in  our  held  experiment  were  partly  due  to  changes  in  innate 
susceptibility.  Feeding  rates  were  also  higher  on  induced  foliage  (Supplementary  Information), 
however,  and  differences  in  variability  between  treatments  in  the  laboratory  were  much  lower 
than  in  the  held,  emphasizing  that  induction  can  affect  transmission  in  nature  through  multiple 
mechanisms. 

Our  held  experiment  thus  showed  that  induced  defenses  can  affect  baculovirus  transmission 
in  nature  in  ways  not  apparent  in  laboratory  experiments.  To  understand  the  consequences 
of  our  experimental  results  for  baculovirus  epidemics  and  insect  outbreaks,  we  used  our  esti¬ 
mates  of  average  infection  risk  D  and  variability  in  infection  risk  C  in  mathematical  models 
(Supplemental  Information).  First,  a  standard  result  from  the  theory  of  epidemics  states  that 
reductions  in  variability  in  infection  risk  can  lead  to  increased  epidemic  intensity12.  Accord¬ 
ingly,  inserting  our  estimates  of  average  infection  risk  v  and  variability  C  into  an  epidemic 
model  shows  that,  rather  than  reducing  the  virus’s  ability  to  cause  population  collapse,  the 
main  effect  of  induced  defenses  is  to  cause  higher  infection  rates  at  high  density  (hg.  ID), 
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in  agreement  with  data  from  epidemics  in  both  gypsy  moths10  and  other  forest-defoliating 
insects9,11. 

To  understand  the  consequences  of  this  effect  for  outbreaks,  we  modified  existing  out¬ 
break  models  to  allow  for  induced  defenses  (Supplementary  Information).  In  current  models, 
outbreaks  are  terminated  by  epidemics  of  a  specialist  pathogen,  such  as  a  baculovirus,  while 
inter-outbreak  populations  are  held  in  check  by  generalist  predators  or  parasitoids4,22.  With 
the  inclusion  of  environmentally  driven  stochasticity,  the  models  can  show  long-period,  large- 
amplitude  cycles  that  recur  at  irregular  intervals,  matching  the  key  features  of  outbreaks  in 
nature.  It  is  often  the  case,  however,  that  models  that  incorporate  the  wrong  mechanism  can 
accurately  reproduce  observational  data23,  and  so  experimental  tests  of  the  models  are  crucial. 
In  particular,  the  models  only  show  realistic  outbreaks  for  intermediate  values  of  variability 
in  infection  risk  C,  but  our  data  show  that,  on  non-induced  foliage,  variability  C  >  1,  which 
is  high  enough  to  guarantee  that  the  host  population  in  the  models  approaches  a  stable,  point 
equilibrium  (fig.  2A),  whereas  on  induced  foliage  C  «  1,  which  is  low  enough  to  produce 
unstable  oscillations  in  the  models  (not  shown,  see  Dwyer,  et  al.  2000).  Neither  a  stable 
equilibrium  nor  unstable  oscillations  is  consistent  with  the  occurrence  of  regular  outbreaks, 
and  our  data  therefore  show  that  existing  models  have  an  important  flaw.  If  we  instead  allow 
induced  plant  defenses  to  drive  changes  in  variability,  however,  the  resulting  models  produce 
stable  cycles  with  a  long  period  and  a  large  amplitude  (fig.  2B),  and  allowing  for  general¬ 
ist  predators  and  stochasticity  produces  cycles  that  are  as  irregular  as  outbreaks  in  nature 
(fig.  2C,  D).  Induced  defenses  may  thus  play  a  crucial  role  in  driving  outbreaks,  by  reducing 
the  stabilizing  effects  of  variability  in  infection  risk  enough  to  allow  cycles  to  occur  (the  evo¬ 
lution  of  host  resistance  can  play  a  similar  role24,  and  we  suspect  that  both  mechanisms  are 
important  in  nature). 

The  theory  of  insect  population  dynamics,  whether  classical25,26  or  recent4,22,  has  long 
focused  on  natural  enemies,  but  our  data  have  demonstrated  that  natural  enemies  may  in¬ 
teract  with  resource  quality  in  ways  that  have  a  significant  effect  on  population  dynamics. 
Meanwhile,  inferences  about  the  effects  of  induced  defenses  that  rely  only  on  laboratory  ex¬ 
periments  have  led  to  conclusions  that  are  contradicted  by  observations  from  the  field27,  but 
basing  models  on  held  experiments  has  instead  allowed  us  to  reconcile  experimental  and  ob- 
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servational  data.  More  concretely,  baculoviruses  are  widely  used  to  control  pest  populations  of 
forest  insects13,  and  previous  work  implied  that  induced  defenses  would  interfere  with  control 
by  reducing  average  infection  risk'.  In  contrast,  our  work  suggests  that  induced  defenses  can 
increase  infection  rates  and  thus  enhance  control  by  reducing  variability  in  infection  risk.  A 
realistic  theory  of  insect  population  dynamics,  which  incorporates  interactions  between  host 
plants  and  natural  enemies,  may  thus  have  an  important  role  to  play  in  forest-pest  manage¬ 
ment. 
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Induction  method 

Treatment 

Pre-treatment  conc’n  (%) 

Post-treatment  conc’n  (%) 

Jasmonic-acid  spray 

Spray 

19.72  ±  0.82 

27.80  ±  0.82 

Control 

20.74  ±  0.98 

18.09  ±  0.62 

Natural  defoliation8 

Defoliation 

23.30  ±  1.0 

27.05  ±  1.4 

Control 

23.36  ±  0.9 

19.54  ±  0.9 

Table  1:  Effects  of  experimental  JA  spray  and  natural  defoliation  on  percent  hydrolyzable 
tannin  concentration  in  red-oak  foliage.  For  the  JA  experiment,  the  interaction  between 
treatment  and  week  after  spraying  was  statistically  significant  at  the  p  <  10~4  level  (see 
Supplemental  Information),  and  natural  defoliation  similarly  had  a  statistically  significant 
effect8.  Pre-treatment  concentrations  were  significantly  different  (experimentals:  t  =  2.70, 
df  =  52,  p  =  0.0093;  controls:  t  =  2.15,  df  =  49,  p  =  0.0365)  because  of  natural  variability, 
but  post-treatment  concentrations  were  virtually  identical  within  a  treatment  (treatments: 
t  =  0.440,  df  =  52,  p  =  0.662;  controls:  t  =  1.311,  df  =  49,  p  =  0.1958  ). 


Figure  Legends 


Figure  1:  Effects  of  induction  on  baculovirus  transmission  and  epidemics.  AIC 
analysis  (see  Supplemental  Information)  showed  that  our  data  pro¬ 
vide  extremely  strong  evidence  that  transmission  on  control  trees 
(A)  is  nonlinear,  while  transmission  on  induced  trees  (B)  is  linear 
(symbols  indicate  data,  and  lines  indicate  best-fit  version  of  equation 
(1)).  (C)  Underlying  distributions  of  transmission  rates  on  control 
branches  (bold  line),  and  induced  branches  (light  line).  For  clarity 
here  we  assume  that  the  distribution  of  transmission  rates  is  log- 
normal,  but  equation  (1)  makes  no  distributional  assumptions28. 
For  induced  branches,  there  is  virtually  no  heterogeneity,  and  so  the 
distribution  is  essentially  a  vertical  line.  (D)  Effects  of  induction  on 
infection  risk  in  full  epidemics,  as  predicted  by  inserting  parameter 
estimates  from  our  analyses  into  an  epidemic  model. 

Figure  2:  Host-pathogen  models  with  and  without  induction,  compared  to 
data  from  gypsy  moth  outbreaks.  The  model  without  induced  de¬ 
fenses  shows  damped  oscillations  (A),  which  are  inconsistent  with 
outbreaks,  whereas  the  model  with  induced  defenses  shows  realistic 
long-period,  large- amplitude  cycles  (B).  These  cycles  are  more  regu¬ 
lar  than  cycles  in  nature  (C) 29 ,  but  allowing  for  generalist  predators4 
produces  cycles  that  are  as  irregular  as  cycles  in  nature  (note  that 
the  data  only  give  information  about  the  timing  of  outbreaks,  not 
the  amplitude,  see  Supplementary  Information  for  more  details). 
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Probability  density  Transmission,  -log(1  -i) 


Figure  1: 
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logl  0(Acres  defoliated)  logl  0(Host  density) 


Figure  2: 
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